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This  report  was  technically  reviewed  by  Dr.  Andrew  Assur  and  Dr.  Donald  E.  Nevel  of  CRREL. 

One  of  the  more  difficult  theoretical  problems  in  ice  engineering  is  the  buckling  of  ice  sheets 
under  horizontal  forces.  A relevant  case  may  be  the  end  of  a pier  or  an  isolated  offshore  structure. 
The  present  report  presents  mathematical  solutions  for  a somewhat  simplified  case  which  eventually, 
by  means  of  computers,  can  be  used  to  give  guidance  for  design. 


CONTENTS 


Page 

Abstract ' 

Preface ii 

Nomenclature v 

Introduction 1 

1 . The  problem 1 

2.  Abstract  of  the  result 4 

Part  I.  Fundamental  solutions 4 

3.  Fuchsian  type  solutions 4 

4.  Contour  integral  solution 5 

5.  Integration  of  the  integral  solution 7 

6.  Fundamental  solutions  for  a = 1 .16 

7.  Fundamental  solutions  for  a = 0 18 

8.  Eigenvalues  for  a = 0 20 

9.  Fundamental  solutions  for  a > 1 21 

Part  II.  Asymptotic  expansions 24 

10.  Asymptotic  expansion  for  0 < a < 1 24 

1 1 . Asymptotic  expansion  for  1 = a = 2 25 

1 2.  Asymptotic  expansion  for  2 ^ fl  = " 25 

Part  III.  Eigenvalues 27 

13.  Range  of  eigenvalues 27 

14.  Eigenvalues  for  the  free-edge  condition 27 

15.  Eigenvalues  for  the  clamped-edge  and  simple-edge  conditions 31 

16.  Deflection 33 

Acknowledgement 36 

Literature  cited 36 

Appendix  A.  Analytical  continuation  at  the  singular  point 37 

Appendix  B.  Tensorial  transformations 39 

Appendix  C.  Comparison  of  the  semi-infinite  plate  buckling  with  the  asymptotic  buckling..  43 


ILLUSTRATIONS 

Figure 

1.  Contour  L on  the  complex  plane  f 6 

2.  Free-edge  condition  roots  x0  in  the  range  1 - v2  = a = 2 28 

3.  Values  of  free-edge  condition  roots  x0  expressed  with  logarithmic  scale  in  the  range 

1 - 1>2  = a = 2 using  a as  a parameter .. .. ^ 28 

4.  Free-edge  condition  roots  xn  (n  = 0, 1 , 2 ...)  as  functions  of  k = \Ja  - 1 in  the  range 

0.01  < a < 00 29 

5.  Behavior  of  the  free-edge  condition  roots  *n  (n  = 0, 1, 2, ...)  in  the  neighborhood  of 

xn  for  k = 0 30 

6.  Clamped-edge  condition  roots  xn  (n  = 0, 1, 2 ...)  as  functions  of  k = y/a  - 1 in  the 

range  0.1  = k = 25 31 

7.  Simple-edge  condition  roots  xn  (n  = 0, 1, 2, ...)  as  functions  of  k = \Ja  - 1 in  the  range 

of  0.1  = k = 25  32 

8.  Behavior  of  the  clamped-edge  condition  roots  x„(n  = 0, 1 , 2, ...)  in  the  neighborhood 

of  xn  = 0 for  k = 0 32 

9.  Behavior  of  the  simple-edge  condition  roots  x0  (n  = 0, 1 , 2, ...)  in  the  neighborhood 

of  xn  =0  for  k = 0 32 


iii 


Page 


Figure  Page 

10.  Normalized  deflection  curves Xq,  x,,  andx2  for  the  free-edge  condition,  when 

a = 17 34 

1 1 . Normalized  deflection  curves  Xq,  x, , and x2  for  the  clamped-edge  condition,  when 

a = 17 34 

12.  Normalized  deflection  curves  Xq,  x1(  andx2  for  the  simple-edge  condition,  when 

a = 17 35 

1 3.  Values  of  H(v,  a)  in  the  range  1 - v2  ~ a = 2 36 


NOMENCLATURE 


a nondimensional  eigenvalue  defined  in  (1.3).  See  (C.1-1)  for  the  physical  meaning. 

D flexural  rigidity . 

F(x)  complex  fundamental  solution  given  in  the  integral  form  by  (4.1 1 ),  in  the  series  form  by 
(5d.5),  and  in  the  asymptotic  form  by  (12.6) . 

fn  ( x ) Fuchsian  type  solutions  defined  by  (3.2),  (3.4),  (3.5),  (3.6). 

k subscript  or  superscript  taking  the  value  1 or  2.  When  k = 1 or  2,  the  upper  or  lower  of 

the  double  sign  ± (or  +)  should  be  taken . 

£0  the  characteristic  length  Defined  by  = (D/y)'/i. 

r radial  distance  from  the  center  of  the  hole.  Used  also  as  an  integration  variable. 

v transformed  unknown. introduced  by  the  contour  integral  (4.1 ). 

w deflection . 

w, , w2  fundamental  solutions . 
x = r/i0  ■ 

x0  the  value  of  x at  the  periphery  of  the  internal  hole. 

xn  the  eigenroots  of  the  determinant  equations  02,  = 0. 

0 = exp(37r/'/4)  = (-  1 + /)/V 2 . 

7 density  of  water.  Used  also  as  Euler’s  constant. 

k = \Ja  - 1 . 

p = VI  - a • 

v Poisson’s  ratio. 
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THE  BUCKLING  PRESSURE  OF  AN  ELASTIC  PLATE  FLOATING 
ON  WATER  AND  STRESSED  UNIFORMLY  ALONG  THE 
PERIPHERY  OF  AN  INTERNAL  HOLE 

Shunsuke  Takagi 


INTRODUCTION 


The  design  of  shore  structures  that  must  withstand  oncoming  ice  sheets  is  one  of  the  important 
themes  of  ice  mechanics.  During  the  last  decade,  field  and  laboratory  measurements  have  been  made 
on  various  kinds  of  structures.10  However,  the  failurS  mechanism  of  floating  ice  sheets  has  not  yet 
been  clarified,  and  the  interpretation  of  the  data  has  not  yet  been  satisfactory. 

One  of  the  failure  mechanisms  of  a floating  ice  plate  is  buckling.  The  only  analytical  study 
presently  available  on  this  subject,  however,  concerns  only  the  buckling  of  beams  on  elastic  founda- 
tions.3 

The  difficulty  of  the  analytical  study  of  this  problem  is  partly  due  to  the  multitudinous  boundary 
geometries  of  the  shore  structures  and  ice  plates.  For  this  study  we  chose  an  ideally  simple  case  and 
carried  out  the  analysis  of  a buckling  failure. 

1.  The  problem 

Let  us  suppose  we  have  a thin  elastic  plate  floating  on  water,  extending  horizontally  to  infinity, 
and  stressed  with  uniform  horizontal  pressure  along  the  periphery  of  an  internal  circular  hole.  We 
are  interested  in  formulating  buckling  pressure  and  deflection  at  the  point  of  failure. 

The  vertical  deflection  w of  an  elastic  plate  that  rests  on  a liquid  and  is  subjected  to  the  middle 
surface  forces  /VXXJ  /V  and  /Vxy  is  governed  by  the  differential  equation 

On4  w + yw  = /V  + 2/V  9 2.*  + N b2-w  (1  1 ) 

v ' xx  dx2  *Y  dxdy  yy  dy2  { ■ 1 


where  D is  the  flexural  rigidity  and  y the  specific  weight  of  the  liquid.5  Let  r be  the  radial  distance 
from  the  center  of  the  hole.  In  our  problem,  the  deformation  is  cylindrically  symmetric  around  the 
center  of  the  hole.  Then  (1.1)  becomes  (see  App.  B) 


d 2 w 
dr 2 


+ N, 


ee 


1 dw\ 
r dr  f 


(1.2) 


where  /0  = [Dly)'/a  is  the  characteristic  length,  and  Nrr  and  Nge  are  the  radial  and  hoop  horizontal 
stresses  in  the  plate. 

Following  the  usual  treatment,5  we  assume  that  the  horizontal  stress  components  Nxx,  Nxy, 
/Vyy  are  in  equilibrium  by  themselves.  Then  they  are  derived  from  a biharmonic  function  0 by 


/V„ 


_ 

by2 


r 


n = 

yy  9x2 


/V  = - d-20 

*v  9x9y 

In  the  general  polar  coordinates  they  become 


/V  =1^+1  d2  0 

rr  r 9r  r2  a02 

r0  9/-  ^ r 90  / 
Naa  =£*-. 


In  our  problem,  0 is  a function  of  z only  and  must  tend  to  zero  when  r becomes  infinite. 
Nrg  = 0,  and  A/rr  and  Ngg  may  be  formulated  as 

^rr  = -«T  'o  r'2 

Nee  = <rf  % r2 

where  a is  a nondimensional  constant  chosen  to  be  positive. 

We  introduce  the  nondimensional  length  x 

.x  = r /31  . 

Thus,  (1.2)  becomes 

dA  tv  + 2 d 3 w _ 1 -a  d 2 tv  + 1 -a  dw  + w _ q 
dx*  x dx2  x2  dx2  x3  dx 

Atx  = <*>,  the  condition 

tv  = 0 


&L  = 0 


must  be  satisfied.  At  x = x0,  where  x0  is  the  value  of  x at  the  periphery  of  the  internal 
consider  three  alternative  conditions: 

(1)  the  clamped-edge  condition: 


hole,  we 


tv  = 0 


o 

dx 


2 
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(2)  the  simple-edge  condition: 


w = 0 

d 2 w + v_  dw_  _ q 
dx2  x dx 

(3)  the  free-edge  condition: 

d2  w + £ dw  _ q 
dx2  x dx 

d I d2  w + 1_  dw\  + a_  dw_  - q 
dx  dx2  x dx  J x2  dx 


(1.8) 


(1.9) 


where  u is  Poisson’s  ratio. 

Eq  (1 ,8)2  and  (1 .9),  are  found  from  M,  = 0.  The  second  equation  of  (1 .9)  is  derived  from 
Kirchhoff’s  condition,5  Q.  + (Mr)(dMrg/dO)  = 0.  The  effect  of  horizontal  stress  is  counted  in 
Qr,  as  explained  below.  In  the  rectangular  coordinates  x,  y,  shears  Qx  and  <?y  are  given  by 
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Q, 


3Afxy 

- + N 

dvv 

+ yV 

a w 

dx 

3 y 

' XX 

dx 

xy 

9 y 

a/wxy 

+ ™yy 

- + N 

dvv 

+ 
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dx 

9 y 

1 xy 

dx 

yy 

3 y 

(1.10) 


We  found  these  equations  by  extending  Hetenyi’s3  one-dimensional  treatment  to  two-dimensional 
treatment.  In  the  polar  coordinates  r,  0,  shears  Qr  and  Qg  are  given  (see  App.  B)  by 


dMt 
d r 


Q,  = -^-L  + ~ 


dMrf) 


ee 


+ 3w  N + 1 9w  m 
+ aT  " r 3 6 "re 


(1.11) 
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+ 2 Mrn  + I + 9*  Nrg  + I NRH. 
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The  constant  a is  the  eigenvalue  to  be  determined  to  satisfy  the  boundary  conditions  at  x = x0. 
The  first  step  for  the  solution  of  this  eigenvalue  problem  is  to  discover,  given  a positive  number  a, 
two  real  functions,  w,  (x)  and  w2  (x),  that  are  the  solutions  of  differential  equation  (1 .5)  and  meet 
the  boundary  conditions  at  x = 00  in  (1 .6).  The  value  of  x0  is  not  yet  specified  at  this  stage,  and 
no  boundary  conditions  are  imposed  in  the  finite  range.  We  call  w-\  and  w2  the  fundamental 
solutions.  We  shall  find  them  later  in  the  following  form: 


+ IW2  ~ f 
1 

The  second  step  is  to  express  the  fundamental  solutions  as  power  series  of  x.  Let  fm  (x)  (m  = 
0,  1,2,  3)  be  the  Fuchsian  type  solutions  of  (1.5)  relative  to  x = 0.  We  shall  find  linear  combina- 
tions 


(r2  + s/r4  - 1 y^l2  + ( r 2 + s/r4  - 1 12  e ^ (r4  - 1 )',/j  dr  (1.1 2) 


wk  (x)  = 


^ km  fm(x ) 


(1.13) 


3 


where/?  = 1 or  2,  by  determining  constants  A km  by  use  of  (1.12).  The  solution  w(x)  is  a linear 
combination  of  the  fundamental  solutions 

w(x)  = A (x)  + B w2(x).  (1.M) 

The  third  step  is  to  solve  the  simultaneous  equations  of  A and  B that  are  found  by  substituting 
(1 .14)  into  the  boundary  conditions  of  (1 .7),  (1 .8),  or  (1 .9).  If  the  algebraic  equation  found  by 
letting  the  determinant  of  the  simultaneous  equations  equal  zero  has  one  or  more  than  one  non- 
negative root,  a root  gives  x0.  Our  problem  is  then  solved. 

2.  Abstract  of  the  result 

The  main  feature  of  the  numerical  result  is  as  follows: 

a.  Under  the  clamped-edge  condition,  eigenvalue  a is  in  the  range  1 < a < °°.  It  does  not  ex- 
tend :o  infinity.  Infinitely  many  eigenvalues  exist  in  the  neighborhood  of  *p  = 0* 

b.  Under  the  simple-edge  condition,  eigenvalue  a is  in  the  range  1 < a = °°.  A single  branch 
extends  to  infinity.  Infinitely  many  eigenvalues  exist  in  the  neighborhood  of  x0  = 0. 

c.  Under  the  free-edge  condition,  eigenvalue  a is  in  the  range  1 - v2  = a = °°,  where  v is 
Poisson’s  ratio  of  the  elastic  plate.  Two  branches  extend  to  infinity.  Note  that  the  range  of  a in 
this  case  extends  to  between  0 and  1 . Infinitely  many  eigenvalues  exist  in  the  neighborhood  of 
*0  = 0. 

d.  Our  numerical  computation  is  complete  in  the  range  0 = a = 626  and  in  the  neighborhood 
of  a = but  incomplete  in  the  intermediate  region  626  < a < We  believe  that  the  range  pre- 
sented in  this  paper  covers  all  the  cases  of  our  interest. 

e.  Buckling  deflections  take  the  shapes  shown  in  Figures  10, 11,  and  12.  In  these  deflections,  curve 
x0  under  the  free-edge  condition  is  observed  frequently  in  laboratory  experiments  and  field  tests. 
Therefore,  we  may  conclude  that  buckling  under  the  free-edge  condition  is  an  important  mechanism 
of  deformation  and  failure. 

PART  I.  FUNDAMENTAL  SOLUTIONS 

The  main  objectives  of  Part  I are  to  find  A km  in  (1.13)  and  to  determine  the  fundamental  solu- 
tions in  series  forms. 

3.  Fuchsian  type  solutions 

Equation  (1 .5)  has  a regular  singularity  at  x = 0.  The  solutions  relative  to  x = 0 are  the 
Fuchsian  type  power  series  of  x.  Their  indicial  numbers  are 

"o  = 0 
"1=2 
v2  = 1 +p 

= 1 -p 

where 

P = VI  - a . (3.1) 

We  express  the  four  solutions  with  the  single  formula  fm  (x) 


4 


M = 2>"m)  X‘'"’ 

n=0 

where  m = 0, 1 , 2,  3,  and 
*}>m)  = 1 

and  the  rest  of  the  coefficients  of,1"*  (n  = 1 ) are  determined  by  the  recurrence  formula 
= l(*'m  + 4n)(|,m  +4n~  2)(*,m  + 4/»  - 1 + 4ft  - 1 + #l)]_1 


The  individual  forms  of  the  solutions  are 


fi).y  (-un 

T(  43”(2»)!  lj/J  + i-(3-|i)jljn  + l-(3  + fi)j 
43"  (2^+1)!  rfn  + j-  (5  - p)|r[  (n  + x (5  + n)\ 


Ht  (3±/i)  rU-(5±p) 


, M-V  <-1)n  , j2v  M'l  I4V  M'l  4 v 

k+1  fa  44"„!  r[  i + 1-  (2  ± m) I rln  + i (3  ± p)  rln  + i 


(5±u)\ 


.x4n*\tu  (36) 


where  k = 1,2,  and  T(  ) is  a gamma  function.  In  (3.6)  we  have  introduced  the  convention  that 
the  upper  or  lower  of  the  double  sign  ± (or  +)  should  be  taken  according  to  k = 1 or  2,  respectively. 
We  use  a subscript  or  superscript  k exclusively  for  this  purpose.  This  convention  is  observed  through 
out  the  paper. 

Differential  equation  (1.5)  has  an  irregular  singularity  at*  = «>.  In  other  words,  the  solution 
relative  to  x - °°  can  be  found  in  the  form 


e'Xx52  Pn 


oo 

where  X satisfies  X4  + 1 =0.  The  series  pn  x'n  in  this  equation  is  asymptotic  and  divergent. 

n=o 

Therefore  this  equation  does  not  provide  any  means  for  determining  in  (1.13). 

4.  Contour  integral  solution 

In  order  to  find  the  fundamental  solutions,  we  must  transform  (1 .5)  by  means  of  the  contour 
integral 


w(x)  = f v(f)  exf  df 


Figure  1.  Contour  L on  the  complex  plane  f.  Points 
A , and  A 2 represent  points  00  0 on  the  respective  branch. 


where  L is  a contour  in  the  complex  plane  of  f that  shall  be  determined  to  let  a solution  of  (1 .5) 
satisfy  the  boundary  condition  (1.6)  at*  = 00 . Following  the  usual  procedure  (ref.  4,  p.  187-188), 
we  arrive  at  the  differential  equation  of  v(f) 


(1  +r»)i£f  + (23+<7)f2^-  +(9  + 3<7)fv  = 0. 

df3  df2 

The  contour  L selected  for  this  solution  is  shown  in  Figure  1 . 

To  find  the  solution  relative  to  f = °°,  let 

t = fir, 

where 

0 = exp(3ff//4). 

Then  (4.2)  becomes 

(,-4_i)  dlil  + io-3  div  + (23  +a)r2%L  + (9  + 3 a)rv  = 0. 
dr 3 dr2  dr 


(4.2) 


(4.3) 


(4.4) 


(4.5) 


This  equation  has  a regular  singular  point  at  r = The  indicial  numbers  Xm  (m  = 0, 1,2)atr  = 
00  are 

Xo  = 3 
X,  = 2 + #i 
X2  = 7-p 

where  p is  given  by  (3.1).  The  solution  vm  ( r ),  corresponding  to  the  indicial  number  Xm,  is 


‘'m  W = 


n=0 


r»  r'x  m'4n 


(4.6) 


where 


6 


(4.7) 


r 


n-1 

*im)  =TT  (4p  + Xm)(4p  + Xm+2)[(4p  + Xm+2-^)(4p  + Xm+2+p))-1. 

p=0 

The  contour  L must  be  such  that  the  point  f = 0 is  a branch  point  of  v(f).  This  condition  is 
satisfied  by  the  series  v,  (r)  and  v2  (r),  as  shown  by  (4.10)  below,  but  not  by  v0  (/•).  The  series 
V]  (r)  and  v2  (r)  can  be  expressed  by  means  of  hypergeometric  series  F ( , ; ; ) as 


vk(r)  = r'2**  F[(2  ±n)/4,  (4  ± /t)/4;  (2  ±n)/2;r~4] . 

The  hypergeometric  series  are  summed  up  by  use  of  the  formula 

F(a,l+a;2j,z)  = 220*1  (1  - z )-»  (1  + VTTF) 

[Ref.  2,  p.  556,  formula  1 5.1.14) . Thus  vk  (r),  where  k = 1 or  2,  simplifies  to 
vk  W = (f4  - 1 )-*  [ (1 12)  (r2  + )]  . 


(4.8) 


(4.9) 


(4.10) 


This  equation  shows  that  f = 0 is  a branch  point.  Formula  (4.9)  can  be  proved  by  showing  that  the 
right-hand  side  satisfies  the  hypergeometric  differential  equation  of  the  left-hand  side  and  that  they 
both  satisfy  the  same  initial  conditions  at/  = 0. 

Suppose  that  vk  (/-)  on  the  branch  A , B of  L in  Figure  1 is  given  by  (4.10).  Then,  vk  (r)  on  the 
other  branch  A 2 B of  L is  given  by 

( r 4 - 1 )-*  [(1/2)  (r2  + ^ . 

Thus  we  find  the  integral  solution 


w-f 


(r2  + vV4  - 1 )'m/2  + (r2  + yV4  - 1 )"/2  J e<Jxr  ( r 4 - 1 )-H  ^ . 


(4.11) 


Then  the  fundamental  solutions  w,  (x)  and  w2  (x)  are  given  by  the  real  and  imaginary  parts  of 

/» 


Ffx)  = w,  (x)  +/'W2  (x). 


(4.12) 


5.  Integration  of  the  integral  solution 

We  shall  integrate  (4.1 1 ) to  a linear  combination  of  f0  (x),  f,  (x),  f2  (x),  and  f3  (x).  The  first 
step  toward  this  goal  is  to  change  the  range  of  integration  in  (4.1 1 ).  We  use  a complex  variable 
/ = 0r,  where  0 is  given  by  (4.4),  to  transform  (4.1 1)  to 

F(x)  = J I (0*2z2  + v/-24  - 1 )‘m/2  + ( 0' 2 z2  + s/-z*  - 1 )"/2  J ezx  0*1  (-/4  - 1 )->4 


(5.1) 


We  change  the  range  of  integration  0 ~ °°  0 to  the  sum  of  the  ranges  p ~ 0 and  0 ~ - °°.  Thus 
(5.1)  becomes 
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ezx  ,-1  0-1  (z4  + }yv,  dz 


(p-*  z 2 + /y/FT 7 r*12  + (/J-2  72  + /vG^+T  )M/2 

(5.2) 


In  the  above  equations,  the  quantities  inside  the  square  roots  are  made  positive  in  the  respective 
ranges  in  order  to  ensure  correct  transformations.  Letting  z = 0r  and  z = - r in  the  ranges  0 ~ 0 
and  0 ~ respectively,  (5.2)  transforms  to  a linear  combination  of  the  normal  range  integrals 


F(x)  = //>,  (x)  +//»2  (x)  - 0 exp(-/irr //4)s,  (x)  - 0 exp(/in//4)  g2  (x) 

where 

(5.3) 

1 

Ak(x)=  y (r2+/v/:T^T)',l/2e<}xr(1-74)'^^ 

0 

(5-4) 

oo 

0k  (*)  = / (^  + V/-4  + 1 )*  "/2  e"xr  (r4  + 1 )-,/j  </r . 

0 

(5.5) 

Expansion  of  exp(0xr)  transforms  (5.4)  to  a power  series  of  x 

hk  (x)  = + 0x  h\k)  + 1-  (0x)2  h ^ + ... , 

(5.6) 

where 

1 

f,M  = J (z-2  +/Vl^T)?M/2  r"  (1  -t4)-*  dr. 

(5.7) 

o 


Integral  (5.5)  transforms,  as  explained  in  the  next  section,  to  an  ascending  power  series  of  x 


01  (x)  = 9^  +g\^x  +0*11,*1+M  +0^*  X2  + ...  (5.8) 

02  M = 0&2>  + ^ii2)  + x +0^2>  *2  + - <5-9) 

Integration  of  (5.7)  will  be  carried  out  later.  Thus  we  find  F(x)  in  the  following  form 

F(x)  = fi0  + S^2>x1'"  + fl,x  + fi|l1)x1+'i  + S2x2 + ...  (5.10) 

where 

B0  = i(htf  > + - /3(e-{**"i/4)  > + e"*i/4  g^)  (5.1 1 ) 

fi,  = i0(h^  + /j(,2))  - (3(e-^MWi/4)  t/S1’  + e^i/4ffi2))  (5.12) 

fi2  = 1 (/7^>  +/?^2>)-/3(e'(»'’ri/4)  <#>  + e*"^4  <#>)  (5.13) 
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(5.14) 


gjk)  = _ pe1 luni/4)  glk)' 

In  this  calculation  we  have  tentatively  assumed  that  0 < a < 1.  Series  are  arranged  in  the  ascend- 
ing order  on  this  assumption.  The  formulas  for  the  outside  of  the  range  0 < a < 1 will  be  derived 
from  the  formulas  in  this  range.  Note  that  the  entries  in  (5.10),  except  fi,  x,  are  the  first  terms  of 
fm  (x),  where  m = 0, 1 , 2,  3.  The  first-order  term,  x,  is  not  contained  in  any  of  fm  [x);  it  is  proved 
that  6,  = 0.  Therefore,  the  entries  in  (5.10)  are  sufficient  to  express  wk  (x)  as  linear  combinations 

of  (*)• 

If  the  original  series  (r)  in  (4.6)  is  not  summed  up  to  a closed  form  (4.1 0),  the  transformation 
of  integral  (5.1)  to  (5.2)  needs  the  analytical  continuation  of  the  integrand.  However,  the  analytical 
continuation  at  the  singular  point  J = 0 is  not  easy.  In  this  particular  case,  this  difficulty  is  com- 
pletely overcome  by  use  of  the  hypergeometric  series.  Further  discussion  is  presented  in  the  Appendix. 

5 a.  Formulas  of  g^ 

To  formulate  g^  we  shall  successively  develop  (5.5)  into  series  of  x.  Integration  of  these 
formulas  will  be  carried  out  in  the  next  section. 

Lettingx  = 0 in  (5.5),  we  find g^: 

oo 

r - / (r2  + \Jr*  + 1 )*^2^  (r*  + I)*1/2  dr.  (5a.1) 

o 

To  find  gff),  subtract^2'  in  (5a.1)from^2  (*)  'n  (5.5), 

g2  (*)  “ 9 f>2*  = f (r2  + \/V4  + 1 W2  (e"rx  - 1 ) (r4  + 1)'1/2  dr. 

0 

Letting  | = rx,  the  right-hand  side  of  this  equation  becomes 

= x1'"  f (£2  + r'2  (e’*  - 1)  (£4  +X4)"1'2  </£.  (5a.2) 

0 

Lettingx  = 0 inside  the  last  integral,  we  find 

g(2)  = 2m/2  f ^-2(e'«-1  )dt.  (5a.3) 

o 

To  find g\2\  we  subtract  the  product  of  x1_M  and  (5a.3)  from  (5a.2) 

92^) -9^-9^  x1-" 

oo 

- «'•*  / 

0 


Letting  £ = rx,  the  right-hand  side  of  this  equation  transforms  to 


(£4  + x4)"1/2  (£2  + £4  + x4)^  - 2“/2  |"*2 


(e-«-1)rf£. 


where 


0,  (r)  = (r2  + n/FTT  )*»  (r4  + 1 )-#  - l*'2  r*4'2 . (5a.5) 

Because  of  the  inequality  1 = (1  -e'u)/w  = 1 -w/2,  the  integrand  of  (5a.4)  is  uniformly  bounded. 
We  can,  therefore,  let  x •*  0 inside  the  integral.  Thus 


ff(2>  - - / <l>)(r)rdr.  (5a.6) 

0 

To  find  , we  substract  the  product  of  x and  (5a.6)  from  (5a.4) 

92(x)  -42)-42)  xUM-g\2)x  = x2  f 0,  (r)r2  -e-X-- ~ tfr. 

1 r x 

Because  of  the  inequality  0 = (1  - u - e*u)/w2  = - 'A,  the  integrand  of  this  integral  is  uniformly 
bounded.  We  can,  therefore,  let  x -*■  0 inside  the  integral.  Thus 


g{2]  = J / 01  (r)r2dr. 


(5a.7) 


To  find  c/P,  we  simply  differentiate^,  (x)  in  (5.5)  with  regard  to x and  let*  = 0 in  the  result. 
Thus  we  find 


oo 

^P  = - (r2  + s/r4  + 1 )~^2  [r/(r4  + 1 )",/4]  dr.  (5a.8) 

0 

To  find  we  use  (5a.1)  and  (5a.8)  to  derive  the  formula 

oo 

gyW-g^-g^x  = f (r2  + n/FTT  )-^2  (e'r*  - 1 + «)  (r4  + 1 )-*  £/r. 
o 

Letting  £ = rx,  the  right-hand  side  of  this  equation  becomes 

oo 

= x1+44  f ($4  + x4)'14  ($2  + VF+x7  )‘m/2  (e-«  - 1 + ?)  (5a.9) 

0 

Letting x = 0 inside  the  integral,  we  find 

oo 

g(J)  = 2-M/2  j |-m-2  (e-t  _ ! + $)  dl  (5a.l0) 

0 

To  find  g^\  we  use  (5a.9)  and  (5a.10)  to  derive  the  formula 
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j»  "■n-.wpwyi-.Hwww  U]  JI.I  mil  mill 


= xi+m  f | ($4 +x4)-^  (?2  + vs4 + x4  )'m/2  - 2-*1/2  r*2 


(e'*-1  +*)</$. 


Letting  £ = rx,  this  transforms  to 


= x2  / <t>2  (r)  r2  e~fX-~  ^ + t/r 

n r*X2 


(5a.11) 


where 


02  W = I'4  + 1 )"*  + nA4  + 1 )^/2  - 2-',/2  r'M'2 . 


Because  the  integrand  of  (5a.1 1)  is  uniformly  bounded,  we  can  let*  -*■  0 inside  the  integral.  Thus 
we  find 


9(2  * = \ / 02  (r)r2dr. 


(5a.12) 


5b.  Evaluation  of 

The  independent  variable  tj  defined  by 


/•2  + y/r4  + 1 = TJ4* 

is  useful  for  the  following  integrations.  This  equation  yields 
r2  = (1  -v)l(2vV2) 

Vi  +Z4  = (1  + tj)/(2tjH) 

/■<*•  = - [(1  + tj)/( 4n3/j)l  tfi?. 

Substituting  these  in  (5a. 1),  we  find 

<#>  = (2V2)-1  = (2V5F)-,r^Ljit)l^ljJi)sin2^. 

Use  of  tj  enables  us  to  integrate  (5a.8)  simply  to 

Use  of  tj  transforms  (5a.6)  to 
t 


,(2)  = _L 


M 


/ 


1 - (1  + T?)  (I-TJ)^2’-1 


(5b. 1 ) 


(5b. 2) 


(5b. 3) 
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By  letting  1+17  = 2 - (1  - 1}),  the  last  integral  divides  into  two  parts, 


= h ~\h 

where 

1 

y 1 = / [1  - (1  -n)"/2]  rfM4)*1  dt] 

0 

which  integrates  to 

/,  = - 4//i  + 2B(n/2, 1 -m/4) 

by  the  partial  integration  of  ij"(4M_1.  The  remaining  integral  /2 

y2  = / in  -r?r/2-  0 -t?)^2’-1]  >?'(m/4)-1 

0 

integrates  to 

y2  = -B(nl2, 1 - p/4). 

Thus  we  find 

0(2>  = At'1  • (5b. 4) 

We  show  the  results  of  (5b.3)  and  (5b.4)  in  a single  formula 

(5b.5) 

Applying  the  partial  integration  of  (5a.3)  integrates;  applying  the  two  times  of  partial 
integration  of  f'*1"2,  (5a. 10)  integrates.  We  show  the  results  in  a single  formula 

gjk)  = ± 2+m/2  [/i(1  ±JU)]-1  r(1  T/1).  (5b.6) 

We  express  <7^  in  (5a.7)  and^1*  in  (5a.12)  in  a single  formula 

(r4  + 1 y'A  (r2  + y/r*TT  )*m/2  _ 2TW2  r-2*M 
The  use  of  rj  transforms  this  to 

4 0 -*?)*<*?. 

Letting  1 + r/  = 2 — (1  — 17),  this  divides  into  two  integrals 
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1 


<7<k>  = _ 

92  8\/2 


/ 1 "(I  +tl)(1  -T ,)-,i(M/2) 


•P  * 1 / 


r2  dr . 


4k)  = (8V2)'1  K}  +(4V2r  K2, 


/c,  = / ji-[(i-Tj)^/2]  V4*4  (i  -nfdn 


which  integrates  to 


by  the  partial  integration  oft?  4 4 . The  remaining  integral  K2 

% 

k2  = f to  -t)*"/2  -(i  -$rW2Mi  r("'5)/4  0 


integrates  to 


Thus  we  find 


,1*1  = [4^1-1  Ml-' 

= - [VST  (i  r (i-^l  r sin  sllfj 


5c.  Integration  of 
Let 

r2  + / Vi  -t4  = f • 


r2  = (1  + f2)/(2?) 

VT7  = /(i  -r)2/(2?) 

rrfr  = [(r2-1)/(4f2)]t/f. 

Substituting  these  in  (5.7),  we  get 

1 n+i  i ji 

/,(k)  = /2-("+1)/2  J f'  2 2 (1+f2)(n-1)/2rff 
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. (1  eosf  (5c- 7) 

because  1^  and  are  real. 

5 d.  Complex  fundamental  solution 
We  get 

s°  - 5^r  r(H^)  <"■" 

by  substituting  (5b. 2)  and  (5c.4)  into  (5.1 1 ).  We  get 

S,  = 0 (5d.2) 

by  substituting  (5b. 5)  and  (5c. 3)  into  (5.12).  We  get 

(5d'31 

by  substituting  (5b.7)  and  (5c. 5)  into  (5.13).  We  get 

fl(k)  = + m±j£  2^/2  exp  (5d.4) 

* pa  4 

by  substituting  (5b.6)  into  (5.14).  Thus  we  find  the  complex  fundamental  solution 

F(x)  = fi0  f0  (x)  + B2  f}  (x)  + fi*1 » f2  (x)  + B™  f j (x).  (5d.5) 

5e.  Fundamental  solutions  for  0 < a < 7 


When  0 < a < 1,  functions  (x)  (/w  = 0,  1,  2,  3)  are  real.  Therefore,  fundamental  solutions 
w,  (x)  and  vv2  (x)  are  found  by  separating  the  coefficients  of  (5d.5)  into  real  and  imaginary  parts. 


Thus 

(x)  = p0f0  (x)  +p,  f2  (x)  +p2f3  (x)  (5e.1 ) 

w 2 (*)  = do  f 2 (x)  + <7i  ^ < x ) + <72  ^3  (*)  (5e.2) 

where 

"73T  r(IH  r(1TJi)  (5t'31 

^ r(Hi)  r(HJ1)  <*■■» 

Pk  = ? mJLtil  2^/2  cos  (5e  5) 

<7k  = + 2^/2  sin  ^ + *0  . (5e.6) 


The  convention  with  regard  tofc  is  observed  in  (5e.5)  and  (5e.6). 
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6.  Fundamental  solutions  fora  = 1 

When  a = 1,  (1.5)  reduces  to 


+ 2.  + * = 0.  (6.1) 

t/x4  * dx 3 


Nevel6  gave  the  Fuchsia n type  solutions  of  this  equation  with  the  notations 


nev0  (x ) = 


nev,  (x)  = 


nev2  ( x ) = 


nel,  (x)  = nev,  (x)  logx  - 


To  find  the  fundamental  solutions,  the  unknown  function  w(x)  must  be  transformed  to  the  un- 
knbwn  function  v(f)  defined  by  the  contour  integral  (4.1).  The  transformed  differential  equation  is 

(1  + f4)  * + 2f3  v = 0 

and  we  find 

„(f)  = (1  + f4)-J4 

Therefore  the  complex  solution  w(x)  for  a = 1 is  given  in  the  integral  form 


efx 

VhF 


dt- 


(6.2) 


The  contour  L is  the  one  shown  in  Figure  1.  To  find  the  fundamental  solutions  in  the  form  of  the 
linear  combination  of  the  nev  functions,  J.  Dieudonne,1  as  explained  in  Nevel,1  expanded  (6.2) 
into  power  series  in  the  neighborhood  of  x = 0,  and  determined  the  first  few  coefficients.  The 
fundamental  solutions  thus  found  are  denoted  here  by  w!^  and  w2 


(*) 


r2(i)  r2(|) 

— nev0  (x)  - —2—  nev,  (x)  + — nev,  (x) 

4x/2 n 0 2y/2  1 2^2?  2 


(6.3) 
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Wj  (x)  = - nev0  (x)  - a-W  nev2  (s)  + nev,  (v)  - (1  - y + log \/7  ) ne^  (x)  (6.4) 


&s/ti 


4-v/tF 


where  y is  Euler’s  constant  0.57721 56. 

We  shall  show  in  the  following  that  w,  (x)  and  w2  (x)  in  (5e.1)  and  (5e.2),  respectively,  give 

lim  w,  (x)  = w?  (x)  + \/2  Wj  (x) 
a -*  1 

lim  w2  (x)  = (x)  -y/2  w2  (x)  . 

a -*  1 

To  show  this,  note  that 

lim  P0  fo  (*)  = (2\/27F  )_1  T2  A-)  nev0  (x) 
a -*  1 ' ' 

lim  <?0  fl  (*)  = (2\Z2tt  T1  r2^  nev2  (x) 
a -*■  1 

lim  fk+,  (x)  = nev1  (x). 
a -*  1 

We  shall  prove,  therefore,  that 

lim  [p,  f2  (x) +p2f3  (x)]  = -V?  (l -7+ logVl"  + ijnev!  (x) +s/2  nel,  (x)  (6.5) 

<7- 1 7 

lim  l?,  f2  (*)  + <7273  Ml  = V?  (l -7+  l°g\/2  - j)  nev,  (x) -y/2  nel,  (x).  (6.6) 

o^1  ' 1 

The  formula 

Pi  f2  M +P2f3  (*) 

in  the  left-hand  side  of  (6.5)  transforms  to 


00 

E *4P+1  lim  Jl-qW  + Cafr,)] 

7HJ  44n  n\  ii-0  M 


p-*0 

by  use  of  (5e.5)  and  (3.6),  where 

rfl  (2±p)l  r|l  (3±p)|  r[l  (5±p)l  „ 

C(n)  = 2+"/2  T(2  + p)  —» p 1 L4  _ LI4 . — 1 

k r|/»  + l-  (2±p)j  r|n  + l (3±p)|  r Jr?  + 1-  (5 ± p)j 

Taking  the  limit,  we  find  that 
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•M+i 


lim  — [- C,  (n)  + C2  (n)]  - 


p -*•  0 


RW  |’—to8*  **  *^r)| 

where  we  have  introduced  the  convention  that  the  summation  ^ disappears  when  n = 0. 

(Z)n  = Z(Z+  1)...(Z  + r»  - 1)  is  Pochhammer’s  symbol  (ref.  2,  p.  246).  This  equation  proves  (6.5). 
The  formula 

q,  f2  (x)  + q2  f3  (x) 

in  the  left-hand  side  of  (6.6)  transforms  to 


oo 

=5 


x4n+1  lim  — [-  Si  (n)  + S2  (n)  1 


44n  n! 


ti~*  0 


by  use  of  (5e.6)  and  (3.6),  where  5k  (n)  is  given  by  replacing  cos  [(3  + p)rr/4]  in  Ck  (n)  with 
sin  [(3  + |i)7r/4] . Taking  the  limit,  we  find  that 

lim  -L  [-5,  (n)  +52  («)]  = 

H-+0  * 


V2 


Hf)n(j)n]  l1-7"'08^  "4  +|jU1+1  +i  +4^-1)j 


where,  by  convention,  the  summation  ^ disappears  when  n - 0.  This  proves  (6.6). 

p=l 

7.  Fundamental  solutions  for  a = 0 
When  a = 0,  (1 .5)  reduces  to 


(JL  +l  <lV 

\dx2  xdx) 


w + w = 0 


as  may  be  derived  by  putting  Nrj  = Oand/Vfl0  = 0 in  (1.2).  We  decompose  this  equation  into  the 
following  two  equations 


T/\ 

V*2  x dx  ) 


w = 0, 


(7.1) 


shown  with  the  double  sign.  The  solutions  of  the  above  two  equations,  satisfying  the  boundary 
conditions  (1.6)  at  x = are 


w = ker  x ± i kei  x 

giving  the  fundamental  solutions  kerx  and  keix.  Thus,  for  a = 0,  we  have 


(7.2) 
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w(x)  = A kerx  + B keix . 


(7.3) 


We  shall  prove  that  w,  (x)  in  (5e.1)  and  w2  (x)  in  (5e.2)  satisfy 

lim  w,  (x)  = >/2  kerx 
a •*  0 

lim  w2  (x)  = \/2  keix. 
a -*•  0 

First  we  note  that 

lim  f0  (x)  = lim  f3  (x)  = berx 
p -*■  1 n -*■  1 

lim  (x)  = lim  f2  (x)  = 4beix. 

p -*•  1 p-*  1 

Denoting 
1 -p  = X 
we  rewrite  (5e.1)  to 

*1-  Jr(2fi)  r(1r)'f»<(1. 4(2 - » 2" ' >ln r(3 - X| t ' f3 

Then,  letting  X -*■  0 

lim  w,  = —2j=  beix+\/2  (log2  - a)  berx  + \/2  lim 
X-»0  2 V : 2 X -*•  0 

By  use  of  f0  in  (3.4)  and  f3  in  (3.6),  we  find  that 


(7.4) 

(7.5) 


- logx  berx  + 


/_  I \n  „4n 


42n  (2n)!2 


which  proves  (7.4). 

We  rewrite  (5e.2)  to 


*2  =ar^2^r<2*'->i!i"T'f3* 


lim  w2  = ^ nj='  berx  -V?  (1  - y + log\/2  ) bei  x + 


X-»0  2 

By  use  of  f,  in  (3.5)  and  f2  in  (3.6),  we  find  that 


2\/2 


/bf , bf7  \ 


lim 
X -*■  0 


\ax  ax/  42n+2[(2n  + i)!] 


1 

p 


which  proves  (7.5). 

8.  Eigenvalues  for  a = 0 

When  a = 0,  no  horizontal  pressure  works  on  the  plate,  and  buckling  should  not  take  place 
under  any  boundary  conditions.  We  shall  prove  below  that  this  is  true  under  the  boundary  condi- 
tions  (1 .7),  (1.8),  and  (1.9). 

The  following  formulas  are  needed  for  the  proof.  Substituting  (7.2)  into  (7.1),  we  find  the 
relations 


and 


ker"x  + x_1  ker'x  + keix  = 0 


ke^'x+x-1  kei'x  - kerx  = 0. 


(8.1) 


(8.2) 


We  shall  prove  that  no  positive  number  x0  can  satisfy  the  clamped-edge  condition  (1 .7).  The 
determinant  of  (1 .7)  is  given  by 


D\  = 


kerx  ker'x 
keix  kei'x 

when  (7.3)  is  used.  Differentiating  £>, , we  find  the  differential  equation 

<?R±_  + *-l  D,  = ker2x  + kei2x. 
dx  1 

Solving  this  under  the  boundary  condition  that  D,  = Oatx  = 00 , we  find 


D,  = -I  / g(ker2  £ + ker2  £)  di, 


which  is  negative  for  any  positive  x,  proving  our  contention. 

We  shall  prove  that  no  positive  number  x0  can  satisfy  the  simple-edge  condition  (1 .8).  The 
determinant  of  (1 .8)  transforms  to 

kerx  (1  -i')x"1  ker'x  + keix 

keix  (1  - y)x_1  kei'x  - kerx 

This  becomes 
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D2  = - 1 j (-(ker2  (•  + kei2  £)  d%  - (ker2x  + kei2x) 

X 


which  is  negative  for  any  positive  x,  proving  our  contention. 

Finally  we  shall  prove  that  no  positive  number  x0  can  satisfy  the  free-edge  condition  (1 .9).  The 
determinant  of  (1 .9)  transforms  to 


(1  - ker'x  + keix  - kei'x 
(1  - r>)x_1  kei'x  - kerx  ker'x 


This  equation  becomes 


D3  = (ker'2x  + kei'2x)  + 1-  £(ker2  i|  + kei2  £)  d£ 

x 

which  is  positive  for  any  positive  x,  proving  our  contention. 

9.  Fundamental  solutions  for  a > 1 

For  a > 1,p  defined  in  (3.1)  becomes  imaginary  and  must  be  replaced  with  p - in,  where 
k = \Ja  - 1 . (9.1) 

To  compute  T(x  + iy ) for  small  y,  we  used  the  formulas 

oo 

|r(x  + /»/r(x)|2  *TT  l1  +y2/(*+")2]~1  (9-2) 

n=0 


and 


ArgT(x  + />)  = y\p(x)  + ^2  f>V(*  + n)  ~ tan'1  b/(x  + n)]l  (9.3) 

n=0  1 J 

(ref.  2,  p.  256).  These  formulas  can  be  proved  by  use  of  Euler’s  formula  for  the  Gamma  function 
(ref.  8,  p.  237). 

Using  these  formulas  for  small  k,  coefficients  of  F(x)  in  (5d.5)  become 


B0  = \ (2n )-’/}  T2  It  U + k2  (4 p + I)'2]-1 

(9.4) 

B2  = l (2tr)'*  T2  (j)  Jf  [1  + «2  (4p  + 3)-2]'1 

(9.5) 

flW  = T (anV'  R exp  |±  ^ ? in  logVI  + /©j 

(9.6) 

where  R and  0 are  defined  by 
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r(2  ?/k)  = R exp(T  /©). 
They  are  given  by 

R = [mca/sinh  (roc)]^ 


9 = k(1  - 7)  + 


[k//»  - tan'1  («//»)]  . 


There  is  a relation 


(9.7) 


fi0  #2  = cosh  (K7r/2) ] — 1 

Equation  (9.4)  for  B0  and  (9.5)  for  B2  were  used  for  k 2 1 1 in  our  computation.  Equation  (9.7) 
for  © was  used  for  k = 5.  Beyond  these  respective  ranges,  the  asymptotic  expansion  of  the  Gamma- 
function  was  used  to  re-evaluate  B0,  B2  and  0. 

Functions  f0  (x)  and  fy  (x)  are  real.  To  decompose  the  complex  function  fk+1  (x)  into  the  real 
and  imaginary  parts,  we  must  transform  the  denominator  of  (3.3)  to 

(An  + 1 ± n)(An  - 1 ± n)(An)(An  ± 7n) 


which  reduces  to 


= 8/j[2n(1<jfl2  + 4 - 5a)  ± in(32n2  - o)l . 
Thus  we  find 


OP 

fk*1  (*)  = (-  1)n  exp (T/pn)  x4"*1*1* 

n=0 


where 


p<0°>  = 1 
f>0  = 0 


p(0)  _ g-2n 


(«!) 


-1 


TT  [p2(4p2  + 1 -|0)2  + (o-1)(ap2 

P=1 


-'/l 


p " fp 

forr»  = 1,  in  which 

tan  fp  = x(4p2  - 8)/|p^4p2  + 1 - ^ j 

where 

“ » < fp  = » • 
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Fundamental  solutions  tv,  (x),  w2  (x),  and  their  derivatives  are  the  real  and  imaginary  parts  of 
f(x)  and  its  derivatives,  respectively.  We  formulated  them  up  to  the  third  derivatives.  The  following 
formulas  apply  when  m = 0, 1,  2,  and  3: 


d m w 


dxm 


1 = ^-(l)  -IT  A f . JL=-1. . V' 

2V25T  'J0  \ (4p  + 1 )2  / 


-1 


t/xm 


_ /?  exp<Kir/4) 

QK 


oo 

(-  1 )"  /i,m) X 4n+ 1 ~m  005(0,  -pn+0<"V 


dm 


f 2 - 


t/xm  yjliia 


R eXEl^KTr/4)  y (_  ! )n  p(m)  x4n+1-m  cos(0  + 0(mh 
(7K  ' ' n ' 2 n ^ n * 

n= 0 

lr(i.  j-v  r 

2TO  fax  (4p  + 3)2  / 


r2(3- 


o-x"1 


R expiKg/4)  ^ (_1)n/>|)m>x4n+1-msin(01 -pn+0<,m)) 


n=0 


+ R_ex£^Knl4i  ^ (- 1 )n /^m)  x4n+1'm  sin(02  + pn  - 


n=0 


where 


Ok  = 


|±Klog-^F(1-7)K?|J  ^-tan-'^j 
p(1)  = p(0)  [(4n  + 1)2+x2],/2 


P<2)  = pjl)  [(4n)2  + k2 J,/,J 
PnV  = p(n2>  [(4/7  - I)2  +K2]’/j 

0(0  = o 

0*,1  * = tan-1  [x/(4n  + 1)] 

42,  = ^,  + f 

0 [,2)  = 0J,”  + tan-1  (k/(4«)] 


for/J  > 1 


(9.8) 


(9.9) 
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for  n > 1. 


0<3>  = 0{,2*  + Jr  - tan"1  k 
0*,31  = 0* ;2)  + tan'1  [k/(4/7  - 1 )] . 

In  our  numerical  computation  these  formulas  were  effective  only  in  the  range  0 = x = 10. 

PART  II.  ASYMPTOTIC  EXPANSIONS 

Values  of  a series  solution  developed  in  Part  I must  overlap  on  a certain  range  of  x with  the  values 
of  an  asymptotic  expansion  determined  corresponding  to  the  respective  series  solution.  The  series 
must  be  used  for  any  x larger  than  the  overlapping  range. 

10.  Asymptotic  expansion  for  0 < a < 1 

We  can  define  the  analytical  contir  .on  of  vk  (r)  in  (4.8)  by  the  analytical  continuation  of  the 
hypergeometric  function  F(  , ; ; ).  Thus  we  can  transform  i/k  (r)  originally  defined  in  the  range 
1 < r < <*>  to  the  one  defined  in  the  neighborhood  of  r = 1 

vk(r)  = „<k>  f(  1 (2  + M),l-(2-p);|;l-r4)  + 

+ rV-1)-*F(l  (2  + n),  i-  (2-p);l;1  -r4)  , (10.1) 

where 

v<k>  = -2v^  r^I(2±f/))  r^I(2±M))  r(±i// 

v<k)  = VtF  r(l  (2  ±p)j  |r(I  (2 ±m))  r(i  (4  ±p)jj _1 . (10.3) 

Double  signs  may  not  appear  in  the  hypergeometric  functions  on  the  right-hand  side  of  (10.1)  be- 
cause of  their  symmetric  properties  with  regard  to  the  first  and  second  parameters. 

Letting r = 1 + t and  developing  the  hypergeometric  functions  on  the  right-hand  side  of  (10.1 ) 
into  a power  series,  we  can  find  a complex-form  asymptotic  expansion  for  0 = a = 1. 


P}  = - (1  + 2</)/4 
P2  = (9  + 200  + 4<72)/96 
S,  = - (3  + a)/6 
B2  = - (3  + o)(5  + c)/1  20. 

Asymptotic  expansions  of  w,  (x)  and  w2  (x)  are  given  by  the  real  and  imaginary  parts  of  (10.4), 
respectively. 

11.  Asymptotic  expansion  for  1 = a = 2 

A form  of  asymptotic  expansion  for  a = 1 is  found  by  letting  pi  = it c in  the  coefficients  of 
v<'>  + v(k2)  (k  = 1,2).  In  this  case,  formulas  (9.2)  and  (9.3)  need  to  be  modified  to  include  the  case 
x = 0.  The  modified  formulas  are 


j r(/»| 2 = y‘2  XT  l1  + Wn'>2y' 

n=1 

(11.1) 

Arg  r(/»  = - 1 sign(y)  -yy  + ^ - tan'1 

(11.2) 

where 

sign(y)  = 1 for  y > 0 

= -1  fory  < 0. 

The  result  of  the  transformation  becomes  extremely  simple: 

^(1)  + „(2)  = cos(K  log  2) 

(11.3) 

and 

+ = k sin|x  + |ogN/2^J. 

(11.4) 

Substituting  these  into  (10.4)  the  complex  form  asymptotic  expansion  for  1 = a can  be  found.  We 
found,  however,  that  this  formula  was  not  effective  for  large  values  of  a.  Therefore,  we  used  this 
formula  only  for  1 f.  a = 2. 

12.  Asymptotic  expansion  for  2 = a = 00 

Letting  p = Ik,  the  integral  solution  (4.11)  transform  to 

j F(x)  = J e0xr  cos  log  ( r 2 + y/r*  - 1 )j  (r4  - 1 )"^  dr. 

(12.1) 

Expanding  the  integrand  in  the  neighborhood  of  r = 1 by  letting  r = 1 + 1,  and  using  the  approxi- 
mations 

log(r2  + V/-4  - 1 ) = 2 \/T  + OW 
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and 


\Jr*  - 1 = 2\ft  + 0(f) 
we  find  the  integral  asymptotic  solution 


F(x ) ~ [ e^x(1+t>  cos(kVT)  . 


To  evaluate  this  integral,  define  the  function 


oo 

Gk  (x)  = l J effxt±iKs/r 

Then  (12.2)  becomes 

F(x)  ~ e^x  [G,(*)+G2(x)]. 

Letting  t = $2,  (1 2.3)  transforms  to 

Gk  w = / exp  (p*  (? ¥ fr)2  + w)d* 


Define  z by 


= -;2. 


The  root  z of  this  equation,  satisfying  the  condition  that  the  real  part  of  z must  approach  positive 
infinity  as  % -*■  °°,  is 


= exp(-|)x^f?exp(5|i.)l 


Use  of  z thus  defined  transforms  (12.3)  to 


exp(5*i/8) 


(1  , \ 5 

^T+T/  I exp(-/2)t/z 


where  we  have  introduced 


ak  = Texp^|  rr/^  1-  xx'%  . 


Transforming  the  contour  of  integration  to  the  sum  of  the  two  contours,  ak  ~ 0 and  0 ~ + °°,  we 
find 

G1(x)  + G2(x)  = sfrfc  exp  (x2/(40x)  + tr//8] . (12.5) 

Thus  we  get 


w-yjt  *i)- 


(12.6) 


This  equation  gives  a little  worse  approximation  in  the  neighborhood  of  a = 1 than  the  one  derived 
in  section  1 1 . 


PART  III.  EIGENVALUES 
13.  Range  of  eigenvalues 

The  clamped-edge  condition  yields  eigenvalues  in  the  range  1 < a < To  show  this,  we 
rewrite  (1.5)  to 


d2  fx  d2  w\ 
dx2  \ dx2  I 


+ xw 


(1  - a) 


A.  Ix^l\ 

dx  [ dx  ) 


(13.1) 


We  integrate  the  product  of  w and  this  equation  between  x0  and  °°;  use  of  the  boundary  condition 
(1 .7)  simplifies  the  successive  partial  integrations.  Thus  we  get 


£7—1 


/ *\{&Y  H%Yd* 


(13.2) 


to  prove  our  contention. 

For  the  simple-edge  condition,  a similar  treatment  yields 


£7-1  = 


-■fr)2];  + / i^r)2  H*|//  (13-3) 


The  right-hand  side  of  this  equation  does  not  seem  to  be  nonnegative;  however,  the  following 
numerical  computation  shows  that  it  is  nonnegative  and  that  1 = a = °°. 

In  the  case  of  the  free-edge  condition,  the  boundary  condition  does  not  simplify  the  partial 
integrations.  The  following  numerical  computation  shows  that  1 -v2  = a = °°.  Eigenvalue  <7 
can  be  less  than  1 only  in  this  case. 

14.  Eigenvalues  for  the  free-edge  condition 
Define  operators 


L = Al.  + — A- 

dx 2 * dx 


M 


d2  .1  d2  1 -£7  d_ 

. ' . — ZL  . • 


£/jf3  X (Jx2  X 2 


dx 


(14.1) 

(14.2) 


Then  the  determinant  D3  for  the  determination  of  A and  B in  (1.14)  is  given  by 

L(w,)  M(w,) 

L(w2)  M(w2 ) 

Root  x0  of  D 3 in  the  range  1 - v2  = a = 2 are  shown  in  Figures  2 and  3. 


°3  = 


(14.3) 
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>8  // 


I0'5  lO'4  10' 3 10' 2 10"' 


*o 

Figure  3.  Values  of  free-edge  condition  roots  xQ  expressed  with  logarithmic 
scale  in  the  range  1 - v2  = a = 2 using  a as  a parameter. 

To  evaluate  small  roots,  we  take  the  first  term  of  each  series  fm  (x)  (m  = 0,  1 , 2,  3)  and  approxi- 
mate F(x)  as  follows: 

F(x)  = B0+B2x2 +bI])x'+ij  + b12)  x'-»  . (14.4) 


With  this  approximation,  the  roots  of  D3  = Oare  given  by  letting  the  real  part  of  L(F)  equal  zero. 
To  prove  this,  rewrite  (14.3)  to  the  complex  form 


m 


i 


Figure  4.  Free-edge  condition  roots  xn  (n  - 0,  1 , 2 ...)  as  functions  of  k = \ja  - 1 in  the 
range  0.01  < a < °°. 


M(F)  = 2 aB2/x. 

Use  of  this  equation  in  (14.5)  proves  our  contention,  because  is  purely  imaginary  for  any  values 
of  a as  shown  by  (5d.3). 

When  0 < a < 1,  we  can  find  in  this  way  that  the  small  roots  x0  are  given  by 


(14.6) 


Because  x0  must  be  nonnegative,  this  equation  shows  that  the  condition  v = n,  i.e.  a = 1 - v2, 
must  be  met.  When  v = \i,  Xg  becomes  zero.  Therefore,  each  curve  in  Figure  2 terminates  on  the 
axis  of  ordinates  at  a = 1 -v*. 

For  small  //,  expansion  in  terms  of  ju  transforms  (14.6)  to 


ln  (tt  x°)  = ? " 7 " i + 3 m2  (- f 1 ^3)  - ^"3  + h ff3)  (1 4J) 

where  f(3)  is  Riemann’s  zeta-function.  Our  numerical  computation  showed  that  (14.7)  gives  a 
close  approximation  over  the  entire  lengths  of  the  curves  in  the  neighborhood  of  a = 1 - 0 in 
Figure  3. 

The  roots  in  the  range  1.01  = a < 626  are  shown  in  Figure  4.  Equations  (9.8)  and  (9.9)  were 
used  to  calculate  Wq,  w2,  and  their  derivatives  in  (14.3). 

The  small  roots  fora  > 1 are  found  by  letting /j  = /7c  in  L[F]  and  taking  the  real  part.  There 
are  infinitely  many  roots  xn  that  can  be  determined  by  solving  the  equation 


tan  log  -j=-  n + - 0 + tan-1  k + tan-1  = exp(x7r/2) 


(14.8) 


If  xn  is  a root  of  the  equation,  xn+1  given  by 
*n+1  = *n  exp(-7r/x) 
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(14.9) 


1 


is  another  root.  Equation  (14.9)  is  true  even  for  an 
extremely  large  or  small  k,  because  (14.8)  holds  for 
any  k. 

More  than  57  roots  are  plotted  in  Figure  4.  Not 
all  of  them  are  connected  with  solid  lines  to  avoid 
confusion.  The  effect  of  v is  significant  only  for  small 
values  of  k in  the  curves  x0  and  x , . 

The  behavior  ofxn  in  the  neighborhood  ofxn  =0 
is  shown  in  Figure  5.  From  (14.8),  we  get 

dx 

lim  = 0.  (14.10) 

n dx 


Therefore,  curves  xn  that  are  approximated  by  (14.8) 
abruptly  drop  in  the  neighborhood  of  xn  = 0,  and 
touch  the  axis  of  ordinates.  This  behavior  is  not  shown 
by  the  curves  x0  and  x, , which  apparently  cannot  yet 
be  approximated  by  (14.8). 

When  a = 1,  we  simplify  (6.3)  and  (6.4)  by  taking 
the  first  terms  of  nev  functions,  substitute  these  approxi- 
mations in  (14.3),  and  find  that  the  small  roots  are 
given  by 

x0  =.  v/2exp[-(1/v) -7  + tt/(2v/2  )).  (14.11) 

= 0 and  M(x  logx)  = 0.  Equation  (14.1 1 ) shows  that 
x0  tends  to  zero,  when  v approaches  zero  as  the  limit. 

The  asymptotic  behavior  of  the  large  roots  can  be  found  by  using  the  asymptotic  complex 
fundamental  solution  F(x)  in  (12.6)  to  compute  L(F)  and  M(F).  Assuming  that  J defined  by 

f = x2/(  4*2)  (14.12) 

is  of  the  ordinary  magnitude  for  large  x,  we  find 

L(F)  ~F(x)(^2_2f  + r2f2) 

M(F)  ~ F(x)  (1  +f2)|exp^  j-f  exp  ^ )J. 

Thus,  we  discover  that  there  are  two  asymptotic  roots 

*h  = x/v^T  (14.13) 


where 


Figure  5.  Behavior  of  the  free-edge  con- 
dition roots  xn  (n  = 0,  7,  2, ...)  in  the 
neighborhood  of  xn  for  k = 0. 


This  calculation  is  fairly  easy  because  M(x) 


Figure  6.  Clamped-edge  condition  roots  xn(n  = 0,  1 , 2 ...)  as  functions  of  k = \/a  - 1 in 
the  range  0.1  S k ^ 25. 


Therefore,  as  a increases  indefinitely,  the  infinitely  many  roots  xn  in  Figure  4 must  converge 
with  either  one  of  the  two  asymptotic  roots  given  by  (14.13).  We  increased  k up  to  25  (i.e.  a = 
626).  As  shown  in  Figure  4,  curve  x0  and  a few  curves  immediately  to  the  left  of  curve  x0  appear 
to  converge  to  the  right-hand  side  asymptotic  branch.  However,  we  did  not  pursue  the  mode  of 
convergence  any  further. 

If  k becomes  larger  than  25,  some  curves  extend  to  the  right  of  x = 10,  which  was  the  limit  of 
the  use  of  the  series  (9.8)  and  (9.9).  When  roots  xn  are  larger  than  10,  we  should  analytically  con- 
tinue series  (9.8)  and  (9.9)  to  a new  range  of  x.  However,  we  did  not  attempt  this,  since  our 
engineering  interest  was  originally  in  small  values  of  a.  We  presently  consider  that  the  mode  of 
convergence  in  the  range  626  < a < °°  is  not  important  to  us. 

15.  Eigenvalues  for  the  clamped-edge  and  simple-edge  conditions 

The  determinants  D , and  D2  to  be  used  in  the  case  of  clamped-edge  and  simple-edge  conditions, 
respectively,  for  the  determination  of  A and  B in  (1.14)  are 


*1 

W1 

w2 

*T 

*-K) 

w2 

L(w2) 

(15.1) 


(15.2) 


The  roots  for  the  clamped-  and  simple-edge  conditions  in  the  range  0.1  = k = 25  are  shown  in 
Figures  6 and  7,  respectively.  Infinitely  many  roots  accumulate,  as  shown  in  Figures  8 and  9,  in  the 
neighborhood  of  x = 0,  when  a is  in  the  range  1 < a = °°  (i.e.  excluding  the  case  a = 1).  The 
clamped-edge  roots  are  independent  of  v,  because  v does  not  appear  in  the  boundary  conditions.  As 
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shown  in  Figure  9,  the  effect  of  v on  the  simple-edge  condition  roots  are  significant  only  for  small 
x0  and  small  x,. 

To  discuss  the  small  roots,  we  express  O,  and  D2  in  complex  forms.  We  find  that  and  D2 
are  equal  to  the  imaginary  parts  of  - FF'  and  - FL(F),  respectively,  where  F is  given  by  (14.4). 
Letting  u = Ik,  we  find  that  the  small  roots  xn  of  0,  and  D2  are  given  by 

tan[x  log(xn/V2)  -8+  (tr/4)  + tan'1*]  = -exp(-xn/2) 


tan[x  log(xn/V2 ) - 0 + (ir/4)  + r]  = exp(-  ktt/2), 

respectively,  where  we  have  introduced  t defined  by 

tanr  = k (1  + k)/(  1 - a + v),  0 = r < r. 

In  both  cases,  therefore,  if  xn  is  a root,  xn+1  given  by 

*n+l  = xn  exp(- jt/k)  (15.3) 

is  another  root.  Because  (14.10)  still  holds  in  these  cases,  curves xn  in  Figures  8 and  9 are  tangent 
to  the  axis  of  ordinates. 

The  roots  for  the  case  a = 1 are  small  and  are  found  by  taking  the  first  terms  of  nev  functions 
in  (6.3)  and  (6.4).  Using  £>,  in  (1 5.1 ) and  D2  in  (1 5.2),  we  find  that  there  is  no  clamped-edge  condi- 
tion root  and  that  there  are  small  simple-edge  condition  roots 

x0  = \/2  exp  [-  (1  lv)  - y - (w/4)  ] . (1 5.4) 

Therefore,  x0  found  under  the  simple-edge  condition  tends  to  zero  when  v approaches  zero  as  the 
limit. 

For  large  k,  use  of  the  asymptotic  complex  fundamental  solution  (12.6)  yields  that  there  is  no 
clamped-edge  asymptotic  root,  and  that  there  is  one  branch  of  simple-edge  asymptotic  roots,  given 
by 

x = k/2  (15.5) 

which  is  shown  on  the  upper-right  corner  of  Figure  7. 

16.  Deflection 

We  calculated  the  deflections  corresponding  to  roots  Xq.Xj,  andx2  by  applying  the  normalization 
w(xn)  = 1 (16.1) 

for  the  free-edge  condition,  and 

max  | w | = 1 (1 6.2) 

for  the  clamped-  and  simple-edge  conditions.  The  cases  for  a - 17  are  shown  in  Figures  10, 11, 
and  12.  All  the  curves  we  calculated  — but  the  cases  for  a = 17  most  conspicuously  — showed  that 
the  suffix  n of  xn  is  the  number  of  bumps  the  curve  xn  must  pass  through  as  x increases  from  x = 
xn  before  curve  xn  assumes  the  same  shape  as  curve  x0. 

Even  when  a is  in  the  range  1 - v2  = a < 1,  curve  x0  in  the  case  of  the  free-edge  condition 
assumes  a shape  similar  to  curve  x0  in  Figure  10.  This  shape  of  deflection  has  often  been  observed 
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Figure  JO.  Normalized  deflection  curves  x0,  and  x-,  for  the  free-edge  condition, 


Figure  1 1.  Normalized  deflection  curves  x0,  x^,  and  x2  for  the  clamped-edge  con- 
dition when  a = 17. 


in  both  laboratory  and  field  when  floating  ice  plates  are  compre-sed.  We  are  now  convinced  that 
buckling  under  the  free-edge  condition  frequently  takes  place. 

However,  the  shapes  shown  in  these  figures  may  be  different  from  the  actual  ones.  First,  the 
maximum  deflection  is  not  necessarily  upward;  in  Figures  10  and  12,  the  normalization  compels  the 
maximum  to  go  upward.  Second,  the  maximum  point  in  the  real  deformation  may  either  go  above 
the  water  and  not  touch  it,  or  go  downward  and  be  underwater.  The  assumption  introduced  for  the 
derivation  of  the  differential  equation  (1 .5)  is  violated  in  these  cases.  The  deflection  shown  in  these 
figures,  therefore,  may  not  be  actual. 


Figure  1 2.  Normalized  deflection  curves  x0,  x , , and  x2  for  the  simple-edge  condi- 
tion, when  a = 17. 


To  estimate  the  size  of  the  fracture  deflection,  we  will  discuss  the  curve  x0  under  the  free-edge 
condition.  The  other  edge  conditions  can  be  similarly  discussed.  The  fracture  takes  place  when  the 
stress  at  xmin  (see  Fig.  10  for  the  definition)  reaches  the  fracture  stress  of.  In  the  general  polar 
coordinates,  stress  components  on,  ogg,  org  are  given  (see  App.  B)  by 


_ Eh  I d2w  . v_  dw  + v_  d2 w\ 

" " 2(1  - v)  \ dr2  r dr  r2  dQ2  ) 

( L + _L  d2  w + r 92w  \ 

V dr  r2  dO2  dr2  I 

h d_  ( ]_  dw  \ 

+ v)  dr  d 0 J 


Eh 


(16.3) 


where  h is  the  thickness  of  the  plate,  E,  Young’s  moduius.  In  our  case  of  axisymmetry,  introducing 
the  nondimensional  length  x defined  by  (1 .4),  the  above  formulas  become 

a = E/z  d2w  + v dw 
" 2(1  -»)tl  dx2  * dx 


uoe 


Eh 


2(1  -l»)8j 


( 1 dw_  + „ d2w\ 
\x  dx  dx2  f 


(16.4) 


ure 


= 0. 


Therefore,  at  point  xmjn,  where  dw/dx  = 0,  we  have 


uoo 


35 


Let  wn(x)  be  the  w(x)  normalized  at*  = x0.  Then 
the  deflection  is  given  by 

w(x)  = Kwn(x)  (16.5) 

where  K shall  be  determined  by  applying  the  condition 
that 


0.6  0.8  1.0  1.2  1.4  1.6  1.8  2.0 


I °rr  I = °f  at  x = xmj» 
where  af  is  the  fracture  strength.  We  find  K: 
K = 28$  Of  (EA)*1  H(v,  a) 


Figure  13.  Values  of  H(v,  a)  in  the  range 
7 -v2  £ a g 2. 


H(y,a)  = (1 


-■Vi 


Values  of  H(i>,  a),  calculated  for  the  case  1 - v2  £ a £ 2,  are  shown  in  Figure  13. 
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APPENDIX  A.  ANALYTICAL  CONTINUATION  AT  THE  SINGULAR  POINT 

We  analytically  continued  vk  (r)  in  (4.8)  across  the  singular  point  z = 0 in  Figure  1 by  use  of 
the  hypergeometric  function  F{  , ; ; ),  as  shown  in  (10.1).  Let  us  now  suppose  that  vk  (r)  is  a 
series  whose  analytical  continuation  is  not  readily  available.  If  the  following  proposition  is  appli- 
cable, we  can  define  the  analytical  continuation: 

Proposition:  Let  0O  (r),  0,  (t),  and  02  (f),  where 

r = r-  1 (A.1) 

be  the  Fuchsian-type  series  solutions  of  the  differential 
equation  (4.5)  in  the  neighborhood  of  r = 1,  i.e., 
x = 0 irt  Figure  A1 . If  a linear  relationship 

vk  (r)  = Ak0o(r) + fik01  (f)  + Ck02(f)  (A.2) 

with  constants  Ak,  Bk,  Ck  holds  true  at  least  at  one 
point  in  the  hatched  crescent  shape  region  in  Figure 
A1 , the  analytical  continuation  of  vk  ( r ) is  defined  in 
the  circle  that  has  center  B and  passes  through  B 1 
and  By 

Proof:  The  hatched  crescent  shape  region  is  the 
common  region  of  convergence  of  the  two  types  of 
series  vk  ( r ) and  0m  ( t ) (m  = 0, 1 , 2).  In  fact,  the  region 
of  convergence  of  vk  (r)  is  outside  the  circle  that  has  a 
Figure  A 1.  The  common  region  of  con-  center  at  0 and  passes  through  B,  B],  fl2,  and  B 3 ; the 

mergence  ofvk  frj  and  0m  (t)  (m  = 0,  7,  region  of  convergence  of  0m  (t)  (m  = 0, 1 , 2)  is  inside 

2)  on  the  complex  plane  ?,  where  0 = the  circle  that  has  a center  at  B and  passes  through  fi, 

(-  1 + i)!\[2.  and  By  If  coefficients  A k,  Bk,  and  Ck  are  non-zero 

at  one  point  in  the  hatched  regions,  these  coefficients 
must  be  constant  at  this  point,  because  all  the  determinants  of  the  simultaneous  equations 

^ = Ak<t>0*  Bk<h^  +Ck02 
v'k  = <4k0o+Bk0]  +Ck^2 
v'k  = A k 0Q  + Bk  <t>'\  + Ck  02 

are  Wronskians.  If  (A.2)  holds  true  at  one  point  in  the  common  region  of  convergence,  (A.2)  must 
be  true  at  any  point  in  the  common  region  of  convergence,  because  both  sides  of  (A.2)  satisfy  the 
same  differential  equation  (4.5)  and  the  same  initial  conditions.  The  proposition  is  thus  proved. 

However,  the  analytical  continuation  defined  by  use  of  the  hypergeometric  function  in  (4.8)  is 
not  the  type  of  (A.2),  as  shown  in  (10.1).  We,  however,  did  not  try  to  learn  how  we  can  define  the 
analytical  continuation  of  vk  (r)  in  (10.1)  when  the  above  proposition  is  not  applicable,  although  we 
must  anticipate  the  need  of  the  analytical  continuation  of  such  cases.  In  this  study  we  were  lucky 
to  have  discovered  that  vk  (r)  is  expressible  with  a hypergeometric  series  that  is  summable  in  a closed 
form. 
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APPENDIX  B.  TENSORIAL  TRANSFORMATIONS 


We  used  dyadic  tensor  expressions  to  transform  tensor  components.  The  dyadic  expressions  can 
be  organized  in  such  a way  that  the  geometric  or  mechanical  meanings  of  the  tensor  notations  are 
clear  at  any  stage  of  the  transformations.  We  shall  explain  the  derivations  of  the  tensor  formulas  we 
introduced  by  this  method  for  use  in  this  paper. 

B1.  Transformations  of  (1.10)  to  (1.11) 

Figure  B1  illustrates  Qx  and  <?y,  the  magnitude  per  unit  length  of  the  shears  acting  on  the  yz- 
and  xz-planes,  respectively.  The  first  step  for  achieving  our  objective  is  to  find  vector  Q that  has 
components  Qx  and  <?y. 

Let  cx  and  cy  be  the  unit  vectors  in  the  x-  and  ^-directions,  respectively.  Let  cn  be  the  unit 
vector  outwardly  normal  in  the  xy-plane  to  the  hypotenuse  AB  (see  Fig.  B2).  Because 

cn-cx  = dy/ds 


and 


cn-cy  = dx/ds 
cn  is  given  by 

cn  ds  = cxdy  + cy  dx.  (i) 

Let  Qs  be  the  shear  defined  per  unit  length  along  the  negative  z direction  on  the  plane  hanging 
downward  from  the  hypotenuse  AB  (Fig.  6? ).  The  shear  Qi  is  giver,  b\ 

Qsds  = Qx  dy  + Qy  dx.  (ii) 

Define  vector  Q by 

Qcx  = <?x  (iii) 

Q-cy  = <?y  (iv) 

which  are  equivalent  to 
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N = ^rr  ur  ur  + Nrt>  K u9  + u0  ur)  + N86  u8  V 


(B.1 1 ) 


respectively. 

We  substitute  (B.9),  (B.10),  and  (B.1 1 ) into  (B.2),  carry  out  the  differentiations  by  use  of  (B.7), 
execute  the  dottings  of  the  unit  vectors  ur  and  ue,  and  identify  the  components  with  (B.8);  thus  we 
find  (1.11). 

B2.  Transformation  of  (1.1)  to  (1.2) 

We  shall  prove  that  in  the  general  polar  coordinates 


92v 
’XX  9*2 


92  w . */  92  w 

'xv  dxdy  yy  dy2 


(B.12) 


transforms  to 


n 9 r2 


2/V  dw  ( 1 dwA  + N ( _L  d2  w + 1 ) (R  1 31 

2N*>  dr  \r  dg2  r dr) 


Expression  (B.12)  transforms  to  an  invariant  form,  N ■ yyw,  where  N and  y are  defined  by 
(B.5)  and  (B.3),  respectively,  Wa  dyadic  operator  (Wilson9),  and  the  double-dotting ab  • • cd 
means  (b  • c)(a  • d).  In  the  polar  coordinates, yyw  becomes 

Carrying  out  the  differentiations  as  indicated  by  (B.7),  the  above  transforms  to 


dr2 


Double-dotting  this  expression  with  N in  (B.1 1),  we  find  (B.13). 

B3.  Derivation  of  (16.3) 

Substituting  eq  (1 .4)  of  Mansfield,5  we  can  transform  the  tensor  expression 
~ ®xcxcx  + °ycycy  + Txy  (cxcx +cycx) 

to  an  invariant  form 

c r~  - (£>/(i  -i>2)|  (yyw+  A Aw) 

where  A is  an  invariant  operator 

A = c i-  -c  i- 
* 9 y > dx  ’ 

which,  in  polar  coordinates,  becomes 
^ ” ur  7ST  _ u»  • 


(B.14) 


(B.15) 
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We  substitute^  in  (B.9)  and  A in  (B.1 5)  in  (B.14),  carry  out  the  differentiations  as  shown  by 
(B.7)  and  find  that  (B.14)  becomes 

<7  = ®rrurur  + O#  (urUtf  + Ue  Ur)  + Oee  Ufl  U« 
with  orr,  aw  and  age  given  by  (16.3),  where  h = z/2. 
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APPENDIX  C.  COMPARISON  OF  THE  SEMI-INFINITE  PLATE  BUCKLING 
WITH  THE  ASYMPTOTIC  BUCKLING 


We  may  be  tempted  to  infer  that  the  semi-infinite  plate  buckling  solution  must  be  close  to  the 
solution  of  an  infinitely  large  ho'e  that  can  be  found,  as  explained  below,  by  use  of  the  asymptotic 
complex  fundamental  solution  F(x)  in  (12.6).  In  the  following  we  shall  show  that  this  inference 
is  correct  in  the  case  of  the  buckling  pressures  determined  for  the  free-edge  condition  - in  fact,  the 
pressures  are  different  by  only  7%  — however,  the  forms  of  deflection  are  different.  For  the  simple- 
edge  condition,  this  inference  cannot  be  justified.  The  clamped-edge  condition  does  not  yield  any 
solutions  for  either  case. 


C.1 . Comparison  of  the  buckling  pressures 

Let -p  be  the  pressure  per  unit  thickness  on  the  internal  hole  in  our  analysis.  By  (1.3),,  we 
have 

P = <n§r'o  (CM) 

where  r0  is  the  radius  of  the  internal  hole. 

By  use  of  D,  x0,  and  x,  this  transforms  to 

p = >d±L  d_ 

*1  ^ 

Therefore,  for  large  x0,  p asymptotically  approaches 


P 


(C.1 -2) 


In  the  case  of  the  free-edge  condition,  using 

fo  = 2-V3 
in  (14.14)  yields 

p ~ 4(2 -s/3)(D/t2)  = 1.0717 (O/Cg).  (C.1-3) 

As  shown  later,  the  free-edge  condition  semi-infinite  plate  buckling  pressure  is  given  by 

p = D/Z2.  * (C.1 -4) 

Therefore,  the  difference  of  the  buckling  pressure  is  only  7%. 

In  the  case  of  the  simple-edge  condition,  use  of  (15.5)  yields 

p = 4 DU2.  (C.1-5) 

However,  the  free-edge  condition  semi-infinite  plate  buckling  problem  yields  only  a solution  that 
continues  to  oscillate  at*  = If  we  may  take  this  as  the  solution,  we  have 
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p = 2 D/il  (C.1-6) 

which  is  half  the  value  ofp  in  (C.1-5). 

In  the  case  of  the  clamped-edge  condition,  there  is  no  asymptotic  root  in  the  case  of  the 
asymptotic  buckling  as  discussed  in  section  1 5,  and  there  is  no  acceptable  solution  in  the  case  of  the 
semi-infinite  plate  buckling  as  discussed  in  the  next  section. 

C.2.  Buckling  of  the  semi-infinite  plate 

Let  they  axis  be  the  edge  of  the  semi-infinite  plate.  We  assume  that  the  uniform  pressure  - p is 
applied  per  unit  thickness  along  the  y axis.  Then,  (1 .1 ) reduces  to 

Dd^+yw  = .pdl^  (C.2-1 ) 

dxA  dx2 

Define  the  nondimensional  length  d by 

* = */S0 

and  the  nondimensional  parameter  b by 
p = 2 byi-l 
where 

C0  = (D/y)V\ 

Then  (C.2-1)  becomes 

dlw.  + 2b  + w = 0.  (C.2-5) 

d\A  dt,2 

In  this  equation 

6=0 

because  we  assume  p = 0.  The  solutions  of  (C.2-5)  are 
wk  = exp  (Xk  £) 

where 

Xk  = - 6 ± yjb2  - 1 

In  the  following,  we  consider  only  the  fundamental  solutions  that  satisfy  condition  (1.6)  at  x = 
When  6 > 1,  letting 

6 = cosh  2p 

we  find  four  solutions 

cosi>$,  sin  I'd,  cos(d /•'),  sin  (d/*'), 

where 


(C.2-2) 

(C.2-3) 

(C.2-4) 
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v = expW. 

These  solutions  do  not  satisfy  the  boundary  condition  (1 .6)  at  x = °°,  and  there  is  no  fundamental 
solution  in  the  case  of  b > 1 . 

When  b = 1 , we  find  two  solutions 


cosx,  sin*. 

Therefore,  there  is  no  fundamental  solution  in  the  case  of  b = 1 . 

When  0 = b < 1 , letting 

b = cos2ju  <c-2-6) 

where  ju  must  be  restricted  to  the  range 

0 < n = ff/4  (C.2-7) 

we  find  the  fundamental  solution,  which  we  write  in  the  complex  form 

W = e-(P+iv)E  (C.2-8) 


where 


0 = sin  /a 


v - cos/i. 

The  real  fundamental  solutions  w,  and  w 2 ate  the  real  and  imaginary  parts  of  W 
W = w,  + /'w2. 

The  general  solution  w is  given  by 
w = /4  w,  + Bwj 

where  A and  B are  arbitrary  constants.  Note  that  0 is  positive  because  of  (C.2-7). 
The  free-edge  condition  for  the  semi-infinite  plate  is  that 


iisy  = o 

d% 2 


d3  w 

d? 


+ 2b 


dw 

dt 


0 


(C.2-9) 


at  £ = 0.  The  latter  is  derived  from  (1.10),,  which  in  this  case  becomes 


-O^L 

dx3 


= 0. 


The  determinant  to  evaluate  the  eigenvalue  is 
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W)  w"'  + 2 bw\ 

D]  = 

vv2  ^2  + 2ow2 

which  we  transform  to  the  complex  form 
W"  W"'+2bW 
" f/*  K'"’  + 2AW'' 

to  facilitate  the  computation,  where  the  upper  bar  signifies  taking  the  conjugate  complex  number. 
Substituting  (C.2-8),  we  find  the  eigenvalue 


H = » „ = 

o z 

Therefore,  (C.2-3)  yields  (C.1-4)  by  use  of  (C.2-4). 
To  determine  the  deflection,  we  must  calculate 


from  which  we  find 


£ = tanffJ- 


Thus,  (C.2-9)  transforms  to 


w = K e~tl2  cos 


(#<*!) 


(C.2-10) 


where  /C  is  an  arbitrary  constant.  The  extrema  of  this  curve  occur  when 


(C.2-1 1 ) 


where  n is  a non-negative  integer. 

The  simple-edge  condition  is  that 

w = 0 

= 0. 

d!-2 

The  determinant  to  evaluate  the  eigenvalue  is 
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d2  = 


w-,  w} 
W2  W2 


which  we  transform  to  the  complex  form 
I W W"  I 


to  facilitate  the  computation.  The  eigenroots  thus  determined  are  outside  the  range  of  (C.2-7),  but 
one  of  them  is  at  the  border 


M = 0. 


The  general  solution  for  this  case 

w = A cos  £ + 5 sin  £ (C.2-12) 

does  not  satisfy  the  boundary  condition  at  £ = If  we  substitute  the  eigenvalue  b = 1 of  this 
unacceptable  solution  into  (C.2-3),  we  find  (C.1-6)  by  use  of  (C.2-4). 

The  clamped-edge  condition  is  that 

w = 0 

w'  = 0. 


The  determinant  to  evaluate  the  eigenvalue  is 
*T 

°3  = 

w2  w2 

which  we  transform  to  the  complex  form 


The  eigenroots  thus  determined  are 
p = it  12  + nit 

where  n is  an  integer.  None  of  them  satisfy  the  condition  (C.2-7). 

C.3.  Asymptotic  deflection 

The  nondimensional  internal  radius  x0  is  evaluated  in  the  case  of  the  free-edge  condition  by  xh 
in  (14.13)  to  be 

x0  = x/v/4(2+VT)  (C.3-1) 

and  in  the  case  of  the  simple-edge  condition  by  x in  (1 5.5)  to  be 
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*0  = */2. 


(C.3-2) 


We  shall  derive  the  asymptotic  buckling  deflections  for  the  free-  and  simple-edge  conditions  by  use 
of  the  asymptotic  complex  fundamental  solution  F(x)  in  (12.6). 

Decomposing  F(  ),  we  find 


w,  = £ cos / 


(C.3-3) 


w2  = R sin  / 


R(x)  = y/itjx  exp[- (x/x/2) -k2/(4V2  x)] 


l(x ) = x/\/2  - k2/(4\/2  x)  + 7t/8. 


The  simple-edge  condition  (1 .9)  yields  the  ratio 


^-Mxn  m(w2)x 


where  L and  M are  operators  defined  by  (14.1 ) and  (14.2),  respectively,  and  the  suffix  x0  signifies 
that  those  attached  with  this  must  be  evaluated  at  x0.  Use  of  F yields 


L(F)  = |/32_2f2+/r2^|F 

M(F)  = |(0-r'rh)3+4fh(/3-/Hfh)jF 


where  fh  is  defined  by  (14.14).  Decomposing  L(F)  and  M(F)  into  the  real  and  imaginary  parts  and 
using  the  relation 


1 - 4x  +x2  = 0 


satisfied  by  x = fh,  we  find  A and  B 


A - K (1  +fh)  :os[/(x0)]  + (1  -fh)  sin[/(x0)]j 
B = K (1  + fh)  sin  [/(x 0 ) ] - (1  -fh)cos[/(x0)]| 


where  K is  arbitrary.  Imposing  the  normalization 


w(x0)  = 1 


K is  determined 


K = \B(x0)(]  +fh)]-1. 


Thus  the  normalized  deformation  wN  (x)  is  given  by 


^(x)  = («(x)//?(x0)|  cos[/(x)-/(x0)l  +3-'^  sin[/(x)-/(x0)l  . 


To  compare  this  with  the  semi-infinite  plate  solution,  let 


■ 


X = x0  + n. 

Assuming  tj  to  be  negligible  as  compared  with  x0,  wN  (x)  transforms  to 
wN  to)  = exP  ’l/v/2  ) (cos  (rj ls/2  ) ? (V5- )"’  sin  (tj/V 2 ) | 


Letting 

u = f/V? 

this  transforms  to 


The  extrema  of  this  curve  occur  when 
= (2/7  + 1)tt 


where  n is  an  integer.  This  curve  is  close  to  the  curve  of  (C.2-10)  but  is  different. 
The  simple-edge  condition  (1.8),  yields  the  ratio 

when  (C.3-2)  and  (C.3-3)  are  used.  Letting 
A = - K sin  (ff/8) 

B = K cos  (jt/8) 

where  K is  arbitrary,  we  find  the  deflection 

w(x)  = KR(x)  sin  [x/\/2  -k2/( 4>/2x)]. 

To  compare  this  with  the  semi-infinite  plate  solution,  let 
x = k/2  + t?. 

Assuming  rj  to  be  negligible  as  compared  with  k/2,  and  letting 
K - V*/( 2rr)  exp  (k/\/2  ) 
we  find 

w(x)  = sin  (\/2  ?}). 

This  is  a case  of  the  general  solution  (C.2-12),  as  may  be  seen  by  letting 


(C.3-4) 


49 
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